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A small mass particle traveling in a curved spacetime is known to trace a background geodesic 
in the lowest order approximation with respect to the particle mass. In this paper, we discuss the 
leading order correction to the equation of motion of the particle, which presumably describes the 
effect of gravitational radiation reaction. 

We derive the equation of motion in two different ways. The first one is an extension of the 
well-known formalism by DeWitt and Brehme developed for deriving the equation of motion of an 
electrically charged particle. Constructing the conserved rank two symmetric tensor, and integrating 
it over the interior of the world tube surrounding the orbit, we derive the equation of motion. 
Although the calculation in this approach is straightforward, it contains less rigorous points. In 
contrast to the electromagnetic case, in which there are two different charges, i.e., the electric charge 
and the mass, the gravitational counterpart has only one charge. This fact prevents us from using 
the same renormalization scheme that was used in the electromagnetic case. In order to overcome 
this difficulty, we put an ansatz in evaluating the integral of the conserved tensor on a three spatial 
volume which defines the momentum of the small particle. 

To make clear the subtlety in the first approach, we then consider the asymptotic matching of two 
different schemes, i.e., the internal scheme in which the small particle is represented by a spherically 
symmetric black hole with tidal perturbations and the external scheme in which the metric is given 
by small perturbations on the given background geometry. The equation of motion is obtained from 
the consistency condition of the matching. We find that in both ways the same equation of motion 
is obtained. The resulting equation of motion is analogous to that derived in the electromagnetic 
case. We discuss implications of this equation of motion. 
PACS number(s): 04.30.Db, 04.25.-g 



I. INTRODUCTION 



The gravitational waves from an inspiralling binary is one of the most promising sources expected to be detected 
by the near-future interferometric gravitational wave detectors such as LIGO, VIRGO and LISA |,|. In order to 
extract the information of the binaries from the last inspiralling stage, we need accurate theoretical templates of the 
gravitational wave forms As an approach, perturbations of a black hole by an orbiting small particle has been 
studied Q. In all the previous works, the rate of change of orbital parameters are assumed to be determined by the 
energy balance. Although this prescription is powerful, the background black hole is required to have a sufficient 
number of Killing vector fields in order to relate the outgoing gravitational waves with the rate of change of orbital 
parameters of the particle. The Kerr black hole is an important example which does not have a sufficient number of 
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Killing vector fields and so rigorous discussions are restricted to the case in which the particle is in a circular or a 
equatorial orbit. Thus the effect of the radiation reaction on the motion of a small mass particle in a general curved 
spacetime is an important target of theoretical investigation. 

On the other hand, the effect of the electromagnetic radiation reaction on the motion of a charged particle in a 
curved spacetime was discussed by DeWitt and Brehme [|[ (Hereafter DB). In the electromagnetic case, the total 
energy momentum tensor composed of the particle and field contributions satisfies the conservation law. Its divergence 
is integrated over the interior of the tube surrounding the particle orbit with infinitesimal length. The part of the 
integration which does not vanish in the small tube radius limit is transformed into the surface integrations over 
the both ends of the tube and over the surface of the tube by using the Gauss' theorem. The integrations over the 
top and bottom of the tube give the definition of the particle momenta at both ends and the difference between 
them represents the change of the momentum during this infinitesimal time interval, which is to be equated with the 
momentum flow given by the integration over the surface of the tube. In this way the equation of motion is obtained. 
As shown below, also in the case of gravitational radiation reactions, we can construct a conserved rank two tensor 
defined on the background spacetime, composed of the matter field and the metric perturbation. However, there is 
an important difference between electromagnetic and gravitational cases. In the electromagnetism, we can consider 
an extended charge distribution which is supported by a certain force other than the electromagnetic field. Thus 
it will be natural to assume that the charge and mass distributions of a point-like particle are not distorted by the 
effect of the radiation reaction. Therefore one may consistently assume that the momentum and the electric current 
of the particle are proportional to the 4- velocity of the particle. Moreover the electromagnetic charge e is not directly 
related to the energy momentum of the particle which is proportional to the mass m. Hence, even if the limit of zero 
particle radius is taken, the divergent self-energy, (oc e^) can be renormalized by shifting the bare mass. In the case 
of the gravitational radiation reaction, it is not possible to consider such an ideal point-like particle because every 
force field universally couples with gravity. Even worse, the role of e in the electromagnetism is also attributed to 
m. Thus the simple renormalization scheme does not make any sense. In order to overcome this difficulty, we put 
an ansatz that the particle momentum, defined in a similar way as in the electromagnetism but without taking the 
limit of the small particle radius, be proportional to the 4-velocity of the particle, which will not be justified within 
this framework. Under this assumption, we obtain an equation of motion in the covariant form with respect to the 
background spacetime which is analogous to that obtained in the case of electromagnetic case. 

In order to develop a more rigorous formalism, we consider an extension of the matched asymptotic expansion 
that has been already developed by many authors (e.g., D'Eath [0 and Thorne and Hartle Q). We assume that 
the internal metric which describes the geometry around the particle is represented by a Schwarzschild black hole of 
small mass in the lowest order approximation. As the particle moves in the curved background, it suffers from the 
tidal distortion. This effect is taken into account by the homogeneous linear perturbations of the black hole. Since 
we know i = and 1 homogeneous perturbations of the black hole are purely indebted to gauge degrees of freedom 
as long as both the mass and angular momentum of the black hole stay constant, we set them to vanish. We also 
assume that the external metric is approximated by the linear perturbations on the background spacetime generated 
by a point-like particle. Then we consider a limited class of coordinate transformation which keeps the meaning of the 
center of the particle unambiguous and match the external metric with the internal one in the matching region where 
the both approximations are valid. Then we find that for an arbitrary orbit the consistent coordinate transformation 
does not exist, and this consistency condition determines the equation of motion, which is no different from that 
obtained in the first approach. 

This paper is organized as follows. In section 2 we explain the covariant expansion method of the tensor Green 
function, which becomes important in both approaches discussed in the succeeding two sections. In section 3 we 
discuss the first approach, i.e., an extension of DeWitt and Brehme's electromagnetic radiation reaction equation of 
motion to the gravitational counter part. Section 4 is devoted to the second approach, i.e., the matched asymptotic 
expansion method to derive the equation of motion of a small black hole. We find both approaches give the same 
equation of motion. In section 5 we discuss implications of the result. Section 6 is devoted to conclusion. 

We suppose the reader is familiar with the concept of 'bi-tensors' and some useful tools developed by DB. Following 
DB, we assign the indices a, (3, 7, S, e, C, 77 for the point on the particle trajectory, z{t), and the indices /i, i^, ^, p, cr for 
the field point, x. For the reader's convenience, the notation and basic formulas are summarized in Appendix A. 



II. METRIC PERTURBATION 

In this section, for the later use, we calculate the metric perturbation, Sgp^i,^ induced by a point-like particle on 
the background metric, g^i,. The background metric is assumed to satisfy the vacuum Einstein equations. Thus, in 
the following calculations, we use the fact that the background Ricci tensor vanishes. As we assume that the particle 
mass, m, is small compared with the background curvature scale, L, we approximate Sg^^ by the linear perturbation 
induced by a point-like particle, h^^, in the whole spacetime region except for the vicinity of the world line of the 
particle. We call the region in which this approximation is valid as the external zone. On the other hand, in the 
vicinity of the world line of the particle, the self gravity of the particle dominates and the metric cannot be described 
by the linear perturbation induced by a point-like particle. We call this region as the internal zone. In this section 
we concentrate on the external zone. The calculation is performed in an analogous manner to the case of the scalar 
and vector perturbations developed by DB. 

Here we consider the linearized Einstein equations. We introduce the trace-reversed metric perturbation. 
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and set the harmonic gauge condition, 



(2.1) 



(2.2) 



where h(x) and tp(x) are the trace of h^''{x) and that of tp^'^{x), respectively, and the semicolon means the covariant 
derivative with respect to the background metric. In this gauge, the linearized Einstein equations become 
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- R''i%ix)iP^P{x) = SwGTt^'ix) + 0(h2), 



(2.3) 



where O(h^) stands for terms quadratic or of higher powers in the metric perturbation. Thus we define the tensor 
Green function G^'^°'^{x, z) which satisfies 



G'"'"'''«;^(x, z) + 2R^'(^%{x)G^P°''^{x, z) = -2g°''^^\x, z)g''^^{x, z) 



5^^\z-x) 



-g 



(2.4) 



of Appendix A and g is the 



where gp"{x,z) is the bi- vector of geodetic parallel displacement defined in Eq. 
determinant of the metric g^^{x). 

First we consider the elementary solution Gt'^"'^ {x, z) which satisfies Eq. ( ^.4|) except at the a{x, z) ^ limit and 
takes the Hadamard form, 



G'i'"'^{x, z) = — ^ ^^'^^ + (:e, z) log \a{x, z)\ + wP'"^^{x, z) 



(2.5) 



where a{x, z) is a half of the square of geodetic interval which is defined in Eq. (1.10) of DB. Its property is summarized 
in Eq. (A6) of Appendix A. The functions u^'^"^ {x , z) , v^'^°'^{x, z) and w'^'^"^ {x , z) are bi-tensors which are regular 
in the cr{x, z) — > limit and uP'^°'^(x, z) satisfies the normalization condition. 



lim u'"'°'^{x, z) = lim 2.g"(''(x, z)g''^'^{x, z). 



(2.6) 



If we put the form ( |2.5| ) into the left hand side of Eq. (2.4), the terms can be classified into three parts. One is 
the terms which contain the factor Xja'^ix.z) manifestly and another is the terms w hich contain log |(T(a;, 2:)|. The 
remaining terms have no singular behavior at the cr(a;,z) — *■ limit. Since the form (2.5) is redundant, we can set 
these three sets to vanish separately: 



{^u^^^^H^, z) - ^^l^uP'^'^^ix, z)) a,^{x, z) = 0, 
yt^'^o^P-^i.^i^x, z) + 2R^'(^%{x)v^P°''^{x, z) = 0, 
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{x,z) + I 2w^''"'3'':(a;,z) 



A-^(x,z) 
A(x,z) 



■'"'"'3 (a;, z) ] a.^x, z) + ^'^^"'^'^.^(a;, z) 



^2R<'{ p{x)u^P°'f\x, z) + (w^'^°'''«;4(a;, z) + 2i?^^%(x)w«''"''(x, z)) a{x, z) = 0, 



(2.7) 
(2.8) 

(2.9) 



where we used the bi-scalar A(x, z) defined in Eq. ( |A^ ) in Appendix A. Equation (2.7) is solved with the normalization 
(2.6) as 



u>''"^^{x,z) = 2-g<^{x,z)-g'''^<^{x,z)^J l\(x,z). 



(2.10) 



The functions v^^""^ [x,z) and w^"^"^ {x, z) are to be determined by solving Eqs. (2.8) and (2.£). The function 
(a;, z) is not needed but the function v'^'^"^{x, z) plays an important role in the following discussion. Although 



it is difficult to find the solution of v^""^^ {x, z) in an arbitrary background spacetime, its explicit form is not required 
for the succeeding discu ssio ns. However it is important to note that v'^'^"'^ {x, z) is uniquely determined. The reason 
is as follows. From Eq. (2^) one finds it satisfies a hyper bolic equation. Hence the problem is if its Cauchy data are 
unique or not. First we note the coincidence limit of Eq. (|2.9D, which gives 



lim w'"'"'3(a;,z) = - lim 2g" ^^{z, x)f p){z, x)R''^''P{x). 



(2.11) 
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Then taking the null limit a{x, z) ^ of Eq. ( ^.9[ ), we obtain the first order differential eq uatio n of v^^^^^x, z) which 
can be solved along the null geodesic. Thus this equation with the boundary condition ( ^.11 ) uniquely determines 



{x, z) on the light cone emanating from z. Therefore the hyperbolic equation (|2.8|) has a unique solution. We 



also mention that v^'^"'^ [x^ z) is divergence free, 

v^"'°'l^,^{x,z) = 0. 

To see this we note the harmonic gauge condition on the Green function requires 

lim v^"'°''^.Jx,z)^Q. 

We also see that the equation for v'^'^"^ ^^{x, z) follows from Eq. (^.^), 



(2.12) 
(2.13) 

(2.14) 



where we have used the fact R'^^ '^^.p — 0, which is proved by contracting the Bianchi identities for the vacuum case. 
Thus we conclude that Eq. ( 2.12| ) holds everywhere. 

The Feynman propagator G^"'^{x,z) can be derived from the elementary solution G'i'^"^ {x , z) by the ie- 
prescription. 
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{x, z) = -i^ f^^r^H^ ^ v^-<^(^[:,, z) \og{a{x., z) + le) + w^^^P{x, z)] 
{2ttY \a[x, zj +ie J 



(2.15) 



The imaginary part of the Feynman propagator Gp'^"^ {x, z) gives the symmetric Green function G'^'^"^ {x, z), from 



which we can obtain the retarded Green function G'^J^'^ {x, z), and the advanced Green function G^'^"'^(x, z) as 



G'"'"^{x,z) = --Ira 



G'P'''^ix,z) 



8tt 



i'""^'^{x, z)d{a{x, z)) - v^'^'^ix, z)e{-a{x, z))] , 



G'tr.fix, z) = 20[S(x), z]G'^''"'5(a;, z), 



^Ret 



Grr(^,^) = 2^?[^,S(x)]G^'^"''(a;,z), 



(2.16) 

(2.17) 
(2.18) 



where S(x) is an arbitrary space-like hypersurface containing a;, and 9\T,{x), z] = 1 — 9[z, '^{x)] is equal to 1 when z 
lies in the past of S(x) and vanishes when z lies in the future. 

Now, using the above obtained Green functions, we compute the trace-reversed metric perturbation Tp'^'^(x) induced 
by a point-like particle whose energy momentum tensor is given by 



T'^'ix) dT v^'{x,T)v^x,T) 



S^^Hx-zir)) 
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(2.19) 



where v^{x, r) = g^a{x, z{t))z°'{t) and i"(r) = dz°' /dT. It is also assumed that the particle mass m is small compared 
with the background curvature scale L. Throughout this paper, we take the unit in which L becomes of 0(1). There 
appear relations whose dimensionality looks wrong but in those cases L is omitted for notational simplicity. 

At this point, we must comment on the reason why we consider the point-like particle. Even in the hnear per- 
turbation, in order to generate a general gravitational field in the external zone, we need to consider a source with 
arbitrary higher multipole moments. However, in the following discussion, we are going to discuss the situation in 
which those higher moments are negligible. Thus we consider this special case here. Later it will become clarified in 
what situation this assumption becomes consistent. 

Now, for the point-like particle, the trace-reversed metric perturbation ip'^'^{x) with the retarded or advanced 
boundary condition is computed as 



V Ret /Ad 



Jx) =±2Gm 



1 



(7(x,z(r)) 

■TRet/Adu{x) 



u^^''o.p{x,z{t))z'-{t)zP{t) 



dr«'^''„0(x,z(r))i"(r)i'3(T) ). 



(2.20) 



=poo 



where TRet/Advi^) is the retarded or advanced time of the particle and is a scalar function which is determined by 



4 



cr (x, z{Tif,,t/Adv)) = 0, (2.21) 
9 (S(a;), z(Tflet)) - ^ (^(T^d.), S(a;)) = 1. (2.22) 

Throughout this paper, we take the convention that the upper sign is assigned for the retarded boundary condition 
and the lower sign is for the advanced one. 

Since the expression (2.20) containing the retarded or advanced time T/jet/A[i-u(2^) is not convenient for the com- 
putations in the succeeding sections, we rewrite it by introducing a new parameterization of the field point x. We 
foliate the spacetime with spacelike 3-surfaces perpendicular to the particle trajectory. More strictly, the 3-surfaces 
are defined as a one-parameter family of r by the relation, (T;ct(a;, z(T))i"(T) = 0. We denote the value of t of the 
3-surface containing the point cc by Ta,. That is 

[a;„(a:,z(T))i"(T)]^^^^ =0. (2.23) 

To distinguish the spatial points on the same 3-surface, we use a-_a{x^z(Tx)) and denote the distance between x and 
^(tk) by 



t(x) := y/2aix,ziT,)). (2.24) 

As we are interested in the region where the linear perturbation is valid, e{x) must be much greater than Gm. 
However, concentrating on the region close to the particle even in the external zone, we can take e{x) is smaller than 
the background curvature scale L, because Gm is assumed to be much smaller than L. Thus we can consider the 
expansion assuming that e^x) is small. 

In the following calculation there appear the higher derivatives of i, such as z and i, where a dot means the 
covariant derivative along the trajectory of the particle. Since we are considering the case in which the deviation of 
particle trajectory from the geodesic is a small correction due to the radiation reaction, we suppose that those higher 
derivatives are at most of 0{e{x)). Thus we define the inverse of the reaction time scale t~^{< e{x)) as the scale of 
z, z and so on. We will see later that « Gm and this assumption is found to be consistent. 

We first consider the time retardation or advancement, (^flet/Adu (2^) j 

'5i?et/A<it,(a;) := Tjiet/Adv{x) - T^- (2.25) 
The time retardation or advancement, 5iif.t/Adv{'^)i is expanded with respect to e{x), 

SRet/Adv{x) = T'^{x)k~'^{x) ^1 T ^'^ix)K~^{x) "z " (tj, )cr;Q (o:, zir^)) 

-^^{^)n-\x)-i'{T.S^^O{e^), 



(2.26) 



as given in E g. (4. 40) of DB ^, where k}{x) = — [o-(x, z(t))]^=t-^. Then the expression for the trace-reversed metric 
perturbation ( 2.20| ) becomes 

^iiiv{x) = ±2Gmgf_,a9i^i3 

{±-^K-^z°'7f ~ 4i("z'') - 2Pa'^z'R^se^°'z'^'^ T 2ei?%^5i'^/ 



=poo 



rfr'i;«^a'/3'(2(r.),^(r'))i"'(r')^^'(^') 

dT'v^^^.p,'''{z{T.,),z{T'))z^'{T')zP'{T')\ +0{e\T;^e). (2.27) 

TOO / 



In the above and in what follows, wc omit the suffix Ret/Adv to -i/j^jy for notational simplicity. The detailed derivation 
of the above formula is provided in Appendix A. Hereafter we also suppress the arguments x, z{tx) and Tx unless 
there arises ambiguity. 

The covariant derivative of the metric perturbation becomes also necessary in section 3. For this purpose, we 
calculate the covariant derivatives of the quantities appearing in Eq. (2.27) such as t^, e{x) and so on. With the aid 



In DB, a-oi{x,z{Tx)) is replaced by ~e{x)nia{Tx)i^i. 
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of the relations given in Appendix A, the derivative of can be calculated by taking the derivative of the both sides 
of the identity (|2.23|) as 



(2.28) 



One must be careful that [f{x,z{Tx))]-^ ^ [f{x,z{T))-fj]r=T^ since Tx is x dependent through Eq. ( 2.23 ). By using 
Eq. (2.28) and the relations given in Appendix A, we obtain 



e 



7"^ - K-^z^zf 



(2.29) 



ir^ix, z(r,))]^^ = l^pguj (5^' - i'i') R'^^s.a'^ + 0{e% 



— ( + iV;/3 z ^ + -R'^p^si^a^'^z' ) + O(e^). 



(2.30) 
(2.31) 

(2.32) 



Then using these results and Eq. ( 2.11 ), after some computations, we obtain the 1st and 2nd derivatives of the metric 
perturbation ip^v{x). 

{ ^K^^z'^z'^a''^ - -z°'zP-p - ± 2z'^°'R'^\'^^z^z^ T 2R°'s^^Pz^z^ 



e e 
2 
e 



-z''"RP\^,z^(j'' + - -z^°'R^\,cz'<z^cJ'''z^ - ^z''z^R-'seCZ^<T'"z 

3e 



T / dr'i;"'^c.'/3'^''(z(r.),z(r'))i" (r')i'^ (r') +0(6\r-^), 

"'TOO / 



(2.33) 



V'm^;«p(2;) = '^Gmg^agvpg^igpS 
( 2 



e 3e e 

-^z°'zPR'',\(j'^a''^ - ^z^°'R(^\r„z^''a'^h'a'^z''i - 4^z°' z^ a'^'^ R^\r,,z'a'^ z"^ 



(2.34) 



We note that among the terms on the right hand side of the above expressions, the terms involving the Riemann 
tensor will not contribute to the equation of motion, as will be shown later. 

Using the above expressions, one can directly check that ■0^'^'?^ + 2R^^'^ pip^f = is satisfied for e 7^ to the order 
in which we are concerned. To the contrary, the harmonic gauge condition is not automatically satisfied. We find 



(2.35) 



Generally speaking, thus obtained metric perturbation satisfies the harmonic gauge condition only when the source 
energy-momentum tensor satisfies the divergence free condition, which, in the present point-like particle case, leads to 
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the conclusion that the source trajectory must be a geodesic. This point is totally different from the electromagnetic 
counterpart, in which the Lorentz gauge condition is related to the electric charge conservation. However this does 
not mean the break down of our formalism. Since z is assumed to be a higher order quantity, the missing of the 
harmonic condition is responsible for the neglection of the higher order perturbations. 

III. DEWITT AND BREHME'S APPROACH 

In this section, we develop our discussion in an analogous way to the electromagnetic counterpart given by DcWitt 
and Brehme. First we derive the 'conserved' rank two symmetric tensor. Integrating its divergence over the interior 
of the world tube surrounding the particle, we derive the equation of motion including the effect of the gravitational 
radiation reaction. 

A. Conservation law 

In the formalism developed by DeWitt and Brehme, crucially important was the conserved energy-momentum 
tensor, which consists of the matter and field contributions. In the case of gravity, the matter energy-momentum 
tensor is divergence free by itself in the sense of the covariant derivative with respect to the total metric, ^^^(x). Thus 
the situation looks different. However, if we choose the background spacctime g^vix), we can construct a quantity 
analogous to the conserved energy-momentum tensor in the electromagnetic case. 

We divide the metric 5^1/ (a;) into the background and the deviation from it, 5gfjtv{x), as 

gtiu{x) = g^lu{x) + Sgi_,„{x). (3.1) 

The background metric is assumed to be a solution of the vacuum Einstein equations. We write the Einstein equations 
with matter source, whose degrees of freedom are represented by ^ symbolically, as 

G^" [g + ^g] = SnGT'^'' [g + 5g, (/.] . (3.2) 
We expand G'^'^[g + 6g] and T'^'^[g + 6g, (p] with respect to Sg^i, as 

C-'ig + 5g] = G^"^^^" + G^^^^^iSg] + 6(2+)"'' [,5g], 

(3.3) 

T'^''[g + <5g,0] =t(o)^"^[0] +T(i+)^''[<5g,0], (3.4) 

where the superscript (n) represents the terms of the nth order in the metric perturbation 5g and the superscript 
{n+) represents the nth or higher order terms. Then the Einstein equations are rewritten as 

T'^'^ldg, 0] := r(°)^^[^] + r(i+)'"'[5g, <t>] - ^G^^+^^''[5g\ = ^G(i)''''[<5g]. (3.5) 

Prom the 0{6g) terms of the contracted Bianchi identities, we find G'^^^'^" -vi^^ = when the background is a solution 
of the vacuum Einstein equations. Here we note again that the semicolon means the covariant derivative with respect 
to the background metric. Thus we obtain the covariant conservation law, 

T'^^.fe</'] = 0, (3.6) 

which is what we needed. 

One can see that if the mass of the particle is small enough and so the metric perturbation Sgfj_^{x) is negligibly 
small, the conserved tensor T^^^lSg, (/)] reduces to T(°)'"^[(/)] which is independent of 6g, which implies T^^^'^''.i,[(j)] = 0. 
Thus, in the lowest order of the metric perturbation, one observes that the point particle moves along the background 
geodesic. 

We shall now specify our consideration to the case in which 6g^u{x) is so small that we can replace 5g^i, by the 
linear perturbation h^v induced by T^^^i^^lcj)]. Then the conserved tensor T'^" can be approximated as 

T^-'idg, 4)] « T^"'[\^, </>] = r(°)''''[(/.] + T(i)^^[h, 0] - -i-G(2)''-[h]. (3.7) 

ottG 

The 1st and 2nd terms on the right hand side of the equation vanish outside the matter distribution, while the 3rd 
term bilinear in h does not vanish anywhere. Therefore it may be interpreted as the gravitational contribution to the 
energy-momentum tensor. The 3rd term is expressed in terms of the trace-reversed metric perturbation as 
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+ igf^-^(^-^ip:-^^^.^ + + |- (3.8) 

The derivation of the above formula is given in the Appendix B. 

B. World Tube Around The Particle 

FoUowing DB, we introduce a further parameterization which distinguishes the points on the same 3-surface pa- 
rameterized by r. It is defined by an impHcit relation between x and the 4 parameters, i.e., t^, e, and O* {i = 1, 2, 3), 

a.a{x, z{tx)) = -enai{T:,)D.\ (3.9) 

3 

where O'fJ* = 1. Here nai{Tx) is a set of orthonormal basis on the 3-hypersurface of at z{tx). It is defined by 

i=l 

n„.(T)i"(T) =0, (3.10) 
<7"'='(z(T))n„,(T)n0j(T) =5„-, (3.11) 

at some and Fermi- Walker transported along the trajectory. In the following, we also use fl" := n"iQ^. The world 
tube of the particle is defined in the same manner as in DB. We consider the 4- volume of the interior of the world 
tube between two 3-surfaces of ri and T2, 

Vtube = {x = x{t, e, ft'); n < t < T2, < e < Ctube} ■ (3.12) 

We define the surfaces of this volume, 

^tube = {x = x{t, Ctube, ^'Y, n < T < T2} , (3.13) 

Scap(T) = {x = x{t, e, < e < etube} ■ (3.14) 

The volume measure is given by dV = —g{x)d'^x. The integral measures on T^tube and Scap(r) have been already 
derived in DB, 

dS^(a;)|s,„,, = {e''K\x)-g^„{x, ^(t))0« + 0{e')) \^^^^^ cfQ dr , (3.15) 

d''^ti{x)\^caAr) = T7r-rz;vr«^^;^'"(a;,2(T))i"(r) 



A(x,^(r)) 



ded^n, (3.16) 



(See Eqs. (4.35) and (4.45) of DB with R°'f^{z) = 0. See also Eq. (1.51) of DB. Note that our measure is with respect 
to the vector, while DB defines it with respect to the vector density.). 

Now integrating the null quantity T^^''.^{x) over Vtube and using the integration by part, we obtain 

0= / dv r(3{z{f),z{Tx))fM^-)^x)'^'"'Ax) 

J Vtube 

= f dY.,{x)rp{z{f),z{Tx))fMrx),x)T>'-{x) 

^tuiie +^cap ("''2 ) — Soop (ti ) 

-/ dV[r0{z{f),z{Tx))fMra.),x)\.^T^^^{x), (3.17) 

where f = (ri + r2)/2. The reason why g°' 0{z(f), z{tx)) ij.{z(tx)tX) is multiplied is to make the integrand to be a 
vector at z{f) and to behave as a scalar with respect to x, so that the integration is done in the covariant manner. 
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C. Gravitational Radiation Damping 



To evaluate the second and third hnes of Eq. ( 3.17 ) , we have to put an assumption on the matter configuration. 
We take the tube radius etube so that the tube is in the external zone. Therefore, the metri c per turbation on the 
tube surface 'Stube can be approximated by that induced by a point-like particle given in Eq. ( ^.20 ). Noting that no 
contribution of matter to T^''[h, 0] exists there, T^''[h, (j>] can be computed as 



An 
i 1 



1 



1 
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^ (^8z^°'V'^\sz^'^n^^ - W.ys^^'z'^h'^n^ + 2g°''^V^Sez"'n^z' - m^°'V''^\sP z^ 



0\z 



'1 1 



where z, i, z are evaluated at t = t^, and we have defined 



^'))i"'(r')i^'(r') 

„f>\\;L,a' Il\-I3' (l\ 



J =poo 



(3.18) 

(3.19) 
(3.20) 



=poo 



Here we have denoted the terms which contain the Riemann tensor by T^'^ ^ which is at most of 0(l/e^). It is not 
necessary to write them down explicitly because they do not affect the equation of motion, as will be shown later. 

Noting that the terms containing odd number of VCs in the integrand vanish when d'^fJ-integration is done, and that 
/ = ^ i^g°'^ + z°'z^), the integration over the surface of the tube T^tube for small 5t := T2 — n is evaluated 



as 



d^,{x)tpiz{f),ziT,))g%{z{T,),x)T'^''[h]{x) 



(3.21) 



The surface integration over the cap Scap requires the knowledge of the metric in the internal zone. We assume 
that the matter configuration and the metric is not perturbed much and is kept spherically symmetric inside the 
radius eo(^ Gm) to sufficient accuracy. In other words, locally the particle behaves as if it were an isolated object. 
Mathematically, we assume that the cap integration for e < eo, which in an unrigorous manner defines the particle 
momentum p"g(T), is proportional to i"(T), 



Scap(r),e<eo 



dS,(a;)r/3(^(T),z(T,))/^(z(T,),x)r^''[h](a;) 

TO(eo,r)5"^(z(f),z(T))i''(T) +Gto^O {■m,T^^ ,etube) ■ 



(3.22) 



We call the matter which satisfies this assumption as an ideal point-like particle. 

Outside the radius eo, the metric is approximated by that induced by a point-like particle, (2.20). Then the cap 
integration for eo < e < etube is evaluated by using Eq. (3.18) to give 



Scap(T),e>eo 



dJ:,ix)g"0iz{f),z{T,))f^{z{T,,),x)T'^''[h]ix) 
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^tube / 



Thus setting 



we obtain 



TO(eo, r) = m(r) + 



(3.23) 



(3.24) 



X 



Scap(T2) — ScapCTl) 



dS,(a;)rM^(^),^(^.))5^(^(^x)>2:)T'"'[h](a:) 



^T^l^) + o } ^"(^) + m(f)i"(f) 

2 Ctufc 



5t + 0((5t2 



(3.25) 



Finally, we consider the remaining volume integral Jy ^ dV[g°' pg^ fj\-^uT^'^ . Since [g°' pg" ~ efi, again we expect 
that the integration for e < eo vanishes for the ideal point-like particle to sufficient accuracy. The integration for 
e > eo is at most of O {Gm^eodT). Thus we may neglect the contribution from the volume integral. 

Putting all results together, we obtain the equation of motion, 



(3.26) 



The leading contributions from the terms proportional to the Riemann tensor which we have neglected are at most 
of the zeroth order in etube (or eo)- As they should be vectors, they must take the form Gw? R^^ ^ pa z^ z'^ or 
Gm^R^i^pa-g^'^z'^ by means of the dimensional argument. Thus all the leading terms disappear and the remaining 
terms become higher order in Ctube (or eo). 

Now we consider the normalization condition gafj{z{T))z" {t)z^ (r) = 0. This gives 



i{r) - Gm4 ^V^^sz^Pz' + ^Vpz^ 



(3.27) 



and it can be integrated as 



m(T) = m|l + (^h^P + 1/^) (^W) 



(3.28) 



where we have defined 



ip(y)f,u{x) := T2Gm / dT'v^^a'(i'{x,z{T'))p' {t')P' {t'), 

J =poo 



(3.29) 



which is the part of the trace-reversed metric perturbation due to the i'^j/q^(x, z) term (i .e, th e so-called tail term) 
in the Green function (see Eq. ( ^.2C| )). Since V'(t))^7 is expected to be of 0(Gm), Eq. ( ^.2^ ) tells us that we can 
consistently replace m(T) by m in Eg. (3.26 ) except for the term, rn{T)z". 
Substituting Eq. ( 3.27 ) into Eq. ( 3.2(: ), we finally obtain 



TOz"(r) = -m.\^z°'z'^Pz^ + g°'^{z)Pz^ - ^g°'\z)z'^P 

-li V(z)i* - lg'^'{z)g^f{z)\ (r) ^(„)^,,,(z(r)), 



(3.30) 



where we have also used the relation Vaf3j{z{T)) = '4'{v)af3;j{z{T)) which follows from Eq. ( 2.11 ). If we impose the 
physical boundary condition with no-inco ming waves from the past null infinity, we should take the upper sign (i.e., 
the retarded boundary condition) of Eq. ( |3.29D for 4'(v)i3j- The meaning of this equation is discussed in section 5. 
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IV. MATCHED ASYMPTOTIC EXPANSION 



In this section we give an alternative derivation of the equation of motion obtained in the previous section in a 
more rigorous way by using the matched asymptotic expansion technique. 



A. Matching Scheme 



To begin with, we state the general concept of the asymptotic matching. We first prepare the metrics in both internal 
and external zones by using different schemes. In the internal zone, we expect that the metric can be described by 
using the black hole perturbation. Namely, we assume that the particle is represented by a Schwarzschild black 
hole in the lowest order of approximation. In the present case, the perturbation is caused by the tidal effects of 
the curvature of the spacetime in which the particle travels. We call this construction of the metric the internal 
scheme. In order to make this scheme valid, the linear extension of the internal zone around the particle must be 
much smaller than the background curvature scale. We use the coordinate {X"} = {T, X*} (a = 0, 1, 2,3; i — 1, 2, 3) 
for the internal scheme and |X|(:= V X'^X'^) is assumed to represent the physical distance scale (in this section, we 
adopt the Minkowskian summation rule on a, 6, • • •, and the Kronecker summation rule on ijj, • • • over the repeated 
indices). In the external zone, as discussed in the previous sections, we expect that the metric is represented by the 
perturbations induced by a particle on a given background spacetime. We call this construction of the metric the 
external scheme. We require that the metrics obtained in both schemes be matched in the matching region of both 
zones, by considering the coordinate transformation between the internal and external zones. Safely, we may assume 
the existence of the matching region as long as Gm ^ L ~ 1. we set the matching radius at \X\ ~ {GmLy^^ by using 
the spatial coordinates of the internal scheme, X^. Then, writing down the metric in the internal scheme, we have 
two independent small parameters and Gm/\X\ in the matching region. The power expansion with respect to 

these two small parameters allows us to consider the matching order by order. 

First we consider the expansion of the internal scheme. Recalling that the perturbation in the internal zone is 
induced by the external curvature which is characterized by the length L, the metric can be expanded in powers of 
\X\/L as 

(0) 1 (1) 1 (2) 

9abiX) = H ab{X) + -H ab{X) + H ab{X) + • • • , (4.1) 

(0) (1) 

where H abiX) represents the unperturbed black hole metric. We expect that H ab{X) will be given by the standard 

(1) , ^ 

linear perturbation of the black hole. Later, we find that the matching condition requires that H ab{X) should 
vanish, which is consistent with the notion that the spacetime curvature is of 0(1/L^). Thus the standard black hole 

perturbation theory applies up to H ab\X). As is known well, the linear perturbations of the Schwarzschild black 
hole Q can be decomposed by using the tensor harmonics, which are classified by the total angular momentum, J . 
The monopole mode (J = 0) corresponds to the mass perturbation. Thus we may set this mode to vanish since it 
is natural to suppose that the change of mass due to the radiation reaction is small. The dipole modes (J — 1) are 
related to the translation and rotation. The translation modes are purely gauge and thus we set them to vanish to 
fix the center of the black hole. As we are considering a non-rotating black hole, we also set the rotational modes 
to vanish. In general, the higher modes contain gauge degrees of freedom as well as the physical ones. However, for 
these higher modes, we do not give any principle to fix the gauge for the moment. Further we expand the metric with 
respect to Gm/\X\ as 

(0) (0) ^ (0) 

H ab{X) = rjab + GmH ab{X) + {Gmf H ab{X) + ■■■ , 

(1) (2) 

1(1) 1(1) Gto(i) (Gm)2 (1) 

1^ ^ (0) L (1) L (2) 

1 (2) ^ ^ I (2) Gm (2) ^ ^ (GmY (2) , , , , 

L'^ L'^ (0) L-^ (1) L'^ (2) 

Note that, from the definitions of the expansion parameters, the component of the metric behaves as 

''H ab ^ \X\^'^-^\ (4.3) 
(") 

In order to write down the external metric in terms of the internal coordinates, we consider the coordinate trans- 
formation from X to (T, X"^) given in the form. 
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a.,^{x,z{T)) = ~T^{T,X). 



(4.4) 



We assume X* = corresponds to the center of the particle, — z"{T), hence J^a = at X' = 0. We suppose that 
the right hand side of Eq. (iA) can be expanded in positive powers of as 



T^{T,X) = UdT)X' + i/„.,(r)X'X-'" + ^U,jk(T)X'X^X'' 



(4.5) 



Ahhough it is possible that more complicated terms such as X'^X^ /\X\ may appear, we simply ignore these kinds 
of terms. We shall find it is consistent within the order of the approximation to which we are going to develop our 
consi derat ion below. Here faii - i„{T) is totally symmetric for ii ■ ■ ■ i„ and is at most of 0(L~("~^^). Using Eqs. (A£) 
and ( All ) in Appendix A, the total derivative of Eq. (4^) gives the important relation, 



r,{z{T),x)dx^^i—{T) 



dT 



{T)X' + - 



2 dT ^ ' 



+ ^Ry^s{ziT))fMT)^iT)f,{T)X^X^+0{\X\')\dT 



r^ + r»,(T)x^' + -r,,fe(T)xJx'= 



+ ^R''f3.,s{z{T))f,iT)f\iT)fkiT)X^X^ + 0{\X\^)\dX' 



(4.6) 



Then, with the aid of Eqs. (4.4) and (4.6), the external metric g^y{x) can be transformed into that written in terms 
of the internal coordinates by the relation. 



gab{X)dX-^dX'' = ~g^,{x)dx^'dxr 
Generally, as the external metric can be expanded by Gm/\X\, we write it as 
~9ab{X) = gab{X) + Gm h ab{X) + {Gmf h ab{X) + • 

(1) (2) 



(4.7) 



(4.8) 



Then Gm h abiX) can be recognized as the linear perturbation on the background gab{X). Further we expand it 

(1) 

with respect to \X\/L as 



(0) 

h , 

(0) 



QabiX) = h ab{X) + - h abiX) + -pT ^ abiX) + 

L (0) 

(0) Gm(i) Gto(2) 



1 (2) 



Gm h ab{X) ^Gm h ab{X) + ^ h ab{X) + ^ h ab{X) + 

(1) (1) L (1) L2 (1) 
{Gmf h ab{X) = (Gm)2 T ab{X) + ? ab{X) + ab{X) + 

(2) (2) L (2) (2) 



(4.9) 



As before, 



(n) 



(4.10) 



(m) (m) 

For brevity, we call h ab oi H ab the (™) component and the matching condition for them as the (™) matching. 

(") (") 

In the matching region (|A:| - (Gmi)i/2)^ the (™) component is of O ((Gm/L)("'+")/2) . The matching condition 
requires that all the corresponding terms in Eqs. ( |4.2| ) and (4.9) should be identical. Then what we have to do is to 
equate the terms of the same power in \X\ to desired accuracy. Thus the condition for the (™) matching is 



E 



(Gm)"' (™') 
^ («') 



^ (Gm)"' ("'') 

-y^^-„ ^ (" ) 



/ (Gm)"+^ vi(m-n) 
1 1^1 



(4.11) 
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B. Geodesic; (g) and (J) Matching 

We begin with the (q) and (J) matchings, which are of 0{{Gm/ L)^) and of 0{{Gm/Ly/^) in the matching region, 
respectively. First we consider the external scheme. In these matchings the external metric is the background itself. 
Here, the necessary order of expansion in \X\ is 0(|X|). Since g i^iv{x)dx^ dx'^ = ga0{z)g" ^i,{z, x)g^ i,{z, x)dx'^dx'^ (see 
Eq. (1.33) of DB), we get 



gfj,u{x)dx^ dx'^ 




dT^ ' dT ^ ' 




{T)U{T) 



dz" 

It 



ft{T)U,{T) + 2r,{T)U,k{T)X'' + O 



\xi 

L2 



dX'dX\ 



o 



\xl 

L2 



dTdX' 



(4.12) 



Comparing the above with Eq. ( [l.S| ) and looking at the dependence on X, one can readily extract out h ab and h at 

(0) (0) 

to the lowest order in Gm/L. 

Next we consider the internal scheme. The (g)-component is trivially given by the flat Minkowski metric. To know 
the (J)-component of the internal scheme, it is better to consider all the (^)-components at the same time. Namely 

(1) . 

we consider the linear perturbation of the black hole H ab- For this purpose, we consider the decomposition of the 
linear perturbation in terms of the tensor harmonics as discussed in Appendix C. Since the timescale associated with 
the perturbation should be of the order of the background curvature scale L, it is much larger than the matching 
radius {GmLy/'^. Therefore the perturbation may be regarded as static. It is known that all the physical static 
perturbations regular on the black hole horizon behave as ^ asymptotically, where J is the angular momentum 

(1) 

eigenvalue. However, in H ab, there exists no term which behaves as ^ \X\"^, [m > 2). Hence, except for gauge 
(1) 

degrees of freedom, H ab contains only J = 0, 1 modes. As mentioned before, we set the perturbation of J = 0, 1 

(1) 

modes to vanish. Thus we conclude that we may set H ab — 0- 
The (o) matching becomes 



1 = 



dz\ , ^ „ / Gm 

dr) + 



o-^(r)/„.(r) 



(rT)-component, 
(Ti)-components, 
(i_;/)-components. 



(4.13) 
(4.14) 
(4.15) 



Equations (4.14) and (4.15) indicate that f°"'{T) are spatial triad basis along the orbit, i.e. 



rk{T)f,{T) = g-^{z{T)) + —{T) — {T) 



O 



Gm 



(4.16) 



The (J) matching becomes 



= 2-(T)-^(m'+(^— U 



o = ^(r)/„„-(T)x^ + r.(r) 



(rT)-component, 
'Gm \X\ 



o = 2/„(,(r)r,),(r)x^+ 



Gm \X\ 



L L 
(ij)-components. 



(ri)-components. 



(4.17) 
(4.18) 
(4.19) 



Then the covariant T-derivative of Eq. ( 4.13 ) and that of Eq. ( 4.14 ) with Eq. ( 4.17 ) result in the background 
geodetic motion. 
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One can see that the internal time coordi nate T becomes a proper time of the orbit from Eq. (4.13) in the lowest 
order in Gm lL. I n the same manner, Eq. (4.17) and the covariant T-derivative of Eq. (4.1£) with (ij)-antisymmetric 
part of Eq. (4.18) give the geodetic transport of the triad /"^(r), 



(4.21) 



Further, from Eqs. ( 4.18| ) and ( 1.19| ), we can see 



\ L L 



(4.22) 



C. Hypersurface Condition; (q) Matching 

We now proceed to the (q) matching, in which the external metric is still given by the background but there 
appear non-trivial perturbations of the internal scheme. Although it is of 0{Gm/ L) in the matching region and 
0((G'm/L)i/2) higher than the remaining (]')-component, we consider it first for the reason which will be clarified 
below. 

(2) 

In order to obtain h at, we expand the external metric in terms of the internal coordinates up to 0(|Xp), i.e., we 

(0) 

have to go one order higher than Eq. (4.12). Then the (q) matching becomes 



^|.. = i?....(.(T))^(T)r.(r)^(r)/^-(r)x^x^ + (^^) 



(rr)-component, (4.23) 
^ H^T^ - \^{T)UMT)X'X'^ + lR.Ms{z{T))^{T)fP,{T)f\{T)fu{T)X^X^ + O 

(ri)-component, (4.24) 
H^j = /„(,(T)r,)fcK7^)^'^' + \Ro.p,MT))r.,{T)f,{T)r,{T)fi{T)X^X' + O ^^"^ 



(ij)-component, (4-25) 

where Eqs. (4.21) and ( 4.22| ) have been used to simplify the expressions. Since there appear terms which describe 
J = 2 perturbations of the internal scheme, the (q) matching does not determine f^ijk without specifying the gauge 
condition of the internal scheme. Here we fix it as follows. 

(2) 

As before, we first set J = 0, 1 modes of H ab equal to zero. Applying the discussion given in Appendix C, we find 

(0) 

the (TT')-component does not contain J — 0, 1 modes. As for the (Ti)-component, it contains no J = mode but 
the J = 1 mode is present. Setting it to vanish to the lowest order in Gm/L, we find 

^(T)/...(T) = 0(— ^j. (4.26) 

Turning to the (ij)-component, we find the term with the Riemann tensor contains neither J — nor 1 mode (see 
Appendix C). Hence the vanishing of J = 0, 1 modes implies 

[MT)rm(T)]j^o,i - O (^^) , (4.27) 

where [• • •] j=o.i means the J = 0, 1 parts of the quantity. 

Let us now turn to the modes higher than J = 1. The (TT')-component contains J = 2 mode, but it just 
determines the physical perturbation of the internal metric. The higher modes are absent from the beginning. The 
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(Ti)-component contains J = 2 and 3 modes. As for the J = 2 mode, it does not involve f^f.. Hence it also gives the 

(2) 

physical perturbation of the internal metric. The J = 3 mode of the (Ti)-component of h ab is given by 

(0) 



1 dz" 

2 dT 



(T)/a<„fe>(T)X^X'= + 



Gm \X\ 

~rT2 



(4.28) 



where < ijk > means to take the symmetric traceless part of the indices (se e Ap pendix C). Since this is a pure gauge 
degree of freedom, we may set it to vanish. Then combining this with Eq. (4.26), we find 



dz°' 

It 



iT).Ujk{T) = O 



Gm 1 



(4.29) 



We do not discuss the higher J modes in the (ij)-component, since it does not give us any information which is 
necessary to derive the equation of motion. 

Then, from Eqs. (|t|), (|^ and (|^, we find 



^(r)a;„ (x(T,X),z(r)) = -^{T)T^{T,X) = O 



L4 



-L 



(4.30) 



to the lowest order in Gm/L. Comparing this with the hypersurface condition of Tx^ Eq. (2.23), one finds that the 
T ^constant hypersurface differs from the =constant hypersur face only by 0{e^) = 0(\Xf). Then following the 
calculations done in section 2 again, one finds they are unaltered even if we replace Eq. (2.23) with 



[a;„(a;,z(T))i"(T)]^^^^=0(e4). 



(4.31) 



The only effect of this replacement is to add 0{i^) to the right hand side of the equivalent formula ( A13) in Appendix 
A. Thus T can be identified to the lowest order in Gm/L. The reason why we have done the (q) matching prior to 
the remaining (5) matching is to establish this equivalence of T and t^- 



D. External Perturbation: (?) Matching 



Now we proceed to the first non-trivial order in Gm/\X\. For this purpose, we must develop the external scheme. 
However, since the time slicing by the internal time coord inate T is now identical to that by Tx in the lowest order in 
Gm/L, we can use the previously obtained formula ( ^.27| ) for the external metric perturbation. 

There remains the (5) matching among the the matchings which becomes of 0((G'm/L)^/^) in the matching region. 
This matching relates the masses of the particle in both schemes. Since this matching is oi 0{{\X\/ L)'^), it is allowed to 
consider as if the background external metric were flat. As is well-known, the linear perturbation induced by a point- 
like particle of mass m in the flat background spacetime is exactly equal to the asymptotic metric of a Schwarzschild 
black hole of mass m in the linear order in m. This fact indicates that the matching gives a consistency condition at 
this order. _ — _ i . 

In order to directly check the consistency, we rewrite Eq. ( p.27| ) in terra s of the intern al c oo rdin ates. From ( 4.10 ) 



we have only to consider the first term of the right hand side of Eq. (2.27). Using Eqs. (4.6), (4.7) and the fact that 



e = ^(T;acr'" = ^jTa{T, X)T°'{T, X), the trace-reverse of Eq. ( 2.27 ) is transformed to give 



(0) / 2 2 \ 

Gm h ab{X)dX''dX' - Gm —dT^ + —dX'dX' , (4.32) 

(1) Vl^l 1^1 / 

which corresponds to the asymptotic form of the Schwarzschild black hole of mass to in the harmonic coordinates. 



E. Radiation Reaction; {{) and (f) Matchings 

There are many components which become of 0{{Gm/ LY) and 0{{Gm/ L)'^/'^) in the matching region. However, 
we are interested in the leading order correction to the equation of motion with respect to Gm/L and we found in 
subsection B that in the lowest order the terms which behave as ^ \Xf' or \X\^ determines the motion of the particle. 
Therefore we consider the (}) and (f) matchings here. 
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In order to perform the ({) and (f) matchings, the calculatfon in obtaining Eq. (4.32) must be done to the hnear 
order in \X\. This can be done in the same way with the aid of the expression of k^x) given in Appendix A. Using 
further the results obtained in subsections B and C, we find the following matching equations: 



Gm (1) 
—r- H TT 

L (1) 



Gm (1) 

T H ij 

L (1) 



Gm (2 



(1) 



dz 



(T) + l| + Gm^{T)^{T)Q^,{T) + O 



(Tr)-component, 



dz" 



Gr. 



Gm (1) 

— T. = '^{T)Um + Gm'^{T)f,{T)Q^p{T) + O ( ( ^ 



L 



dT 

(rj)-component, 

{r.{T)UAT) - + Gm/".(T)/,(T)e«^(r) + O 



(ij)-component, 



-jtHtt^2—{T) 



dT 
+Gm< 



dT 



(T)X' 



'^{T)^{T)f\{T)Q^p,{T)X^ 



-U,{T)r,ki{T)X'X^X''X' 



%X 



(4.33) 



(4.34) 



(4.35) 



Gm (2) 
(1) 



R^^^,{z{T)) — {T)f,{T) — {T)f,{T)X^X^ \+0 



^Gm 



L 



(rr)-component, 

^(T)/„,,(r)x^- + r,(r) 



£I^{T)X^ 



where 



+Gm< 



{T)ff{T)f]{T)e^0,{T)X^ - 2R^p,s{z{T)) — {T)f,{T) — {T)f,{T)X^ 



dT 
2 



dz° 



-Ro.f^-,s{z{T)) — {T)f,{T)f\{T)f\{T)X^X' +0 



(Tj)-component, 



dT 
Gm 



L 



(4.36) 



(4.37) 



GmQaf}{T) h(.,)ai3{z{T)), 
GmQap^iT) := h(^^)„p.^^{z{T)), 



(4.38) 



with 



(4.39) 



Note that h(^y^^^(x) is the metric perturbation due to Vf^i^afiix, z) in the Green function. The (ij)-component of the 
(1) matching is not presented here since it will not be used in the following discussion. 

(1) . — -, . 1 

According to the discussion given in subsection B we have H ab = 0. Thus the right hand sides of Eqs. (4.33), (4.34) 

(1) 

I . (2) 

and (4.35) must vanish. As for H ab, we cannot set them equal to zero. However, again from the condition that there 
(1) 



are no J = 0, 1 modes, the corresponding parts in the right hand sides of Eqs. (4.36) and (4.37) must vanish if they 
are extracted out. Following the discussion given in Appendix C, one finds the terms involving the Riemann tensor 
contains no J = 0, 1 mode. Furthermore, the J = 0, 1 mode of the term involving f"jki in Eq. (4.36) is proportional 
to 



(4.40) 
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which van ishes at th e lowest order in Gm/L by Eq. (4.27). Hence the J = 1 modes of the remaining terms in 
Eqs. (4.36) and (4.37) are extracted out to give 



= ^^iT)^iT) + Gm'^{T)^{T)f\{T)Q^p,{T) + O 

(Tr)-component, 
= f.[^{T)^{T) + GmQ„p,{T) — {T)f^mP,m + O \ 

(rj)-component. 



Gm 



L ^ 
L , 



(4.41) 



(4.42) 



The J = mode of Eq. ( 4.3(S| ) is shown to vanish, while that of Eq. (4.37) exists but it does not contain useful 
information. The matching condition for this mode just gives the equation which determines {dz" /dT)faii to the first 
order in Gm/L. 

The covariant T-derivative of Eq. ( 4.33 ) and that of Eq. ( [4.34 ) with Eq. ( 4.41 ) give the equation of motion with 
leading correction due to the radiation reaction, 



n fl^a Gm 

^^(r) = -— (e"^,(T) + Q\p{T) - Q0,%T)) ^{T)'^{T) + O 



V L 



L 



(4.43) 



The leading correctio n to t he evolutio n of t he 'triad' basis, /"i(r), are also obtained from Eq. (4.41) and the covariant 
T-derivatives of Eq. ( |4.35| ) with Eq. ( |4.42| ) as 



D 



rm = (e"^,(T) + e%Mr) - e^.^iT)) f,iT)^iT) + o 



fGm 



L 



(4.44) 



Since the internal time coordinate T is not properly normalized in the external metric, we define the proper time, 
T = t(T), such that (dz/dr)^ = —1. It is easy to see that we should choose 



dT 



G 



dz" , , dz<^ 



-^ = \^ ^eMr)-j-ir) — {r) +0 — 



dT 



dT 



dT 



Gr 



L 



(4.45) 



Since the second term on the right hand side of this equation is proportional to the small perturbation induced by 
the particle, it is guaranteed to stay small even after a long time interval compared with the reaction time scale 
Tr. = {{Gm/L)-^L). Then Eq. ( p^ ) becomes 



L 



(4.46) 



It is easy to see that this equation is identical to that obtained at the end of the previous section. 

Also, the triad basis are not properly normalized in the external metric. Thus we define, e"i(T), such that 
em(r)e"j(r) = (5^ and e",(T) = (5»j +Sjj)/"j - Gm{dz°' /dT){dz^ /dT)p ^Qp^, where s^j is of 0{Gm/L) and the last 
term is added so as to satisfy the orthonormal condition, eai{T){dz°' / dT){T) = (see Eq. (4.34)). We find 



Gm 



Q.p{T)f\{T)f,{T)+0 



(4.47) 



and again this is guaranteed to stay small. Then the evolution of the normalized triad e°'i{T) becomes 



D 



Gm /dz°' dzl^ ^ dz^ 



dT 



^e\{T) = I :^--^e\^-i- + g'^^z)^-^e\ - g-(z)e'^.^ ) (r) e,,.(r) + O 



a.5 ( 



dz^ 



dT dT 



dT 



dT 



dT 



(4.48) 
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V. IMPLICATIONS 



In this section, we first consider the physical meaning of the equation of motion obtained in the preceding two 
sections. Since the equation of motion we have obtained contains an unknown function Vf^^afiix, z), we need to give 
a method to exphcitly determine the particle trajectory. Here we consider a couple of possibilities to calculate it in 
the case of a specific background, such as Kerr spacetime. 

In order to make the meaning of the equation of motion manifest, we divide the perturbed metric in the external 
scheme into two pieces as 

hf,^{x) = h(^)^^{x) + /i(„)^^(a;), (5.1) 



where /i(u)pi/(a;) is the part due to the tail term (defined in Eq. ( 4.39| )) and /i(u)/ji/(2;) is the part due to the u^j^nafi term 



in the Green function (corresponding to the first term on the right hand side of Eq. (2.2C). The singular behavior of 
the perturbed metric in the coincidence limit is totally due to h(^y^)^^{x). Thus, we introduce the regularized perturbed 
spacetime defined by 

9(v)t,,.{x) ■■= 9t,u{x) + ft.(„)^^(a;), (5.2) 
which has no sing ular behavior any more. Then we find the equation of motion ( [4.43| ) and the evolution equation 



of the triad basis ( [4.44[ ) coincide with the geodesic equation and the geodetic parallel transport, respectively, in this 
regularized perturbed spacetime, g(v)^v To see this let us consider the parallel transport of a vector A" along a 
geodesic x°' = z"(r) in this spacetime. It is given by 



to the linear order in h(^y^^^ where 



5Y{v)°'lil ■= 2 (^(f)"/3;7 + h{v)°'j;l3 - ^(t))/37'") ■ (5.4) 



Then one recovers Eqs. ( 4.43| ) and (4.44) by identifying f with T and replacing A" with dz'^/dT or /",;. This fact is 



the main result of this paper. 

We note that the present equation of motion is analogous to that in the electromagnetic case, except that the 
instantaneous reaction force which is proportional to higher derivatives of the particle velocity is absent in the present 
case. This is because the particle traces the geodesic in the lowest order of approximation. If an external force field 
exists, the assumption of the geodesic motion in the lowest order breaks down and furthermore the contribution of the 
external force field to the energy momentum tensor must be taken into account. Since this fact makes the problem 
too complicated, it is beyond the scope of the present paper. 

Now let us consider the way how to construct g(v)fiiy As stated in the beginning of this section, in order to evaluate 
the particle trajectory explicitly, a practical scheme to calculate g(v)fiiy must be developed. Unfortunately, we do 
not have any scheme which can be satisfactorily applied when calculating the particle trajectory with the effect of 
the radiation reaction, even on a specific background spacetime such as a Kerr black hole. Here we just give a few 
primitive discussions on this matter. For definiteness, we focus on the case of the retarded boundary condition. 

Basically, there seems to be two approaches for calculating g{v)tii> (or equivalently h(^^)^^). The first one is to 
calculate /i(t,)^,y directly. The second one is to calculate h^i, = ft.(ti)^iiy + ^(i))/^!/ and subtract from it. In the 

following, we discuss only the first approach. As for the second approach, we have nothing to mention here, but this 
direction of research may be fruitful. 

By definition, evaluated on the particle trajectory is independent of the past history of the particle.p*] 

Therefore if we consider the metric defined by 

hit^\x) = (^<5/<5.- - \g,u{x)g^%x)^ dr'G^-*^(x, z(r'))i"(T')i^(r'), (5.5) 

for any finite Ar (> 0), it will not contain when it is evaluated on the particle trajectory. The difference 

between h^j^'^^ and /i(t,)^i,y comes from the integral over the small interval. 



** There is a possibility that the future light cone emanating from z crosses the particle trajectory. Since inclusion of this 
possibility makes the problem too complicated, we do not consider it here. 
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Gm / dT'vp,^p{x,z{T'))z''{T')z^{T'). (5.6) 



-At 



(Ar) 



Since Vpa-a(3{x, z) is a regular function, this integral will be negligible for a sufficiently small At. Thus liniAr^o 
wiU give 

In the case of the electromagnetic (vector) Green function, a calculation along the above strategy was performed 
by DeWitt and DeWitt |^ by assuming the background gravitational field is weak so that its metric is given by the 
small perturbation on the Minkowski metric, 

Qt.i' = rip^ + h^^l . (5.7) 

DeWitt and DeWitt calculated the relevant part of the Green function perturbatively to the first order in by 
using the Minkowski Green function. An analogous calculation seems possible in the case of gravity to evaluate 

though it seems difficult to develop such calculations to higher orders in h'j^l. 

As an alternative approach, we give below a speculation of a new method in which we do not have to assume the 
deviation from the Minkowski spacetime to be small. The standard perturbation formalism of a Kerr black hole gives 
us the retarded Green function schematically in the form,[^ 

GXp(^^^) ~ E / d^rjT^e— (*^-*^)/.iV™'-(r.,0„^,)/i^'J'-(r.,0„^.), (5.8) 



where hji"^''^{r^9,ip) represents the spatial mode functions with appropriate boundary conditions and Wj^rn.uj is the 
Wronskian. Here we stress the following point. Provided that x is inside the future light cone emanating from 
z, the integrand becomes exponentially small as the imaginary part of uj goes to —oo. Thus the contour of the uj 
integration can be closed in the lower half complex plane oicu. Since the integrand has poles at the quasi- normal mode 
frequencies, at which Wj.m,u} vanishes, the w-integral can be transformed into the summation over the contributions 
from the residues of the poles. We denote this expression of the retarded Green function by 

Cy^(-'^)--2- E (^^Ir') ' e-^-(*^-*^)/^i~(r.,0.,^.)/.i~(r.,0.,^.), (5.9) 

j,m,n ^ <^—<^n 

where w„ are the quasi-normal mode frequencies. On the other hand, when x coincides with z, the replacement of 
the ^-integration into the summation over the quasi-normal modes is not allowed. Therefore lim G'^^*^^'^\x, z) is 

not equal to lim Gf?^*fl(x,z). In particular, we expect that G^^l}9'^\x, z) is regular in the coincidence limit x z, 

since {dWj^m,uj/du;)uj=uj„ is known to behave as l/nl |T^. Thus it may give the regularized retarded Green function 
which we want to obtain. ^ ^ 

Before closing this section, we give comments on some proposals to the equation of motion including the effect of 
the gravitational radiation reaction. One is the use of the radiative Green function (a half of the difference between 
the retarded and advanced Green functions) in the case of a Kerr background [O As seen easily by using the results 
obtained in the previous section, the use of the radiative Green function instead of the retarded one results in the 
replacement of tp{v)fj,u{x) by i'^'^'l^^ix), which is defined by 

/ + 00 
dT'vp^^,^. (x, z(T'))i"' ir')z^' (r'). (5.10) 
-oo 

Gal'tsov has proved that the back reaction force computed using the radiative Green function correctly gives the loss 
of the energy and the z-component of the angular momentum of the particle in quasi-periodic orbits. However, we do 
not think that this fact indicates the correctness of the prescription because those constants of motion are special ones 
which reflect the existence of the corresponding Killing vector field. For example, it is still uncertain if the radiative 
Green function is useful in evaluating radiation reaction effect on the Carter constant. We would like to come back 
to this point in future publication. 

Recently Ori has argued also in the case of a Kerr background that the partial wave decomposition of the retarded 
Green function with respect to the spheroidal harmonics as schematically shown above would result in the automatic 
removal of the divergence when averaged over several orbital periods vM, thus making it possible to derive the 



^^Note that this expression is only schematic, just for the present explanation. The correct expression can be fond in 
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radiation reaction to the Carter constant. We do not think his argument works. Let us consider expansion of the 
Newtonian potential of a particle off the origin of the coordinates by using the spherical harmonics. If one calculate 
the self-force on the particle, one finds the contribution from each partial wave is finite. However, the net self-force is 
intrinsically ill-defined and if one sums up all the contributions from the different modes, it is easy to see that thus 
obtained self force diverges. Thus the finiteness of the partial wave contributions does not imply that of their sum. 



In this paper, we have derived the equation of motion of a particle on a given background including the effect of the 
gravitational radiation reaction, i.e., with corrections of order (Gm) where m is the mass of the particle. Although we 
use the terminology 'radiation reaction' here, it may not be appropriate in a strict sense because the equation of motion 
we have derived may well contain something more than just the radiation reaction. In fact, in the electromagnetic 
case, the existence of the effect which can be termed as 'the induced polarization force on the background spacetime' 
is reported by many authors The existence and the physical meaning of it in the gravitational case is left for 

future study. 

We have derived the equation of motion in two different ways. First we have considered an extension of the 
electromagnetic counter part developed by DeWitt and Brehme Due to the nature of gravitational interaction, 
there appears an ambiguity concerning the renormalization of mass of the particle, which cannot be resolved within 
this approach. We have then derived the equation of motion by setting an ansatz which seems physically reasonable 
without justification. In order to overcome this problem, we have developed a method based on the matched asymptotic 
expansion, assuming the local geometry around the particle is described by a spherically symmetric black hole plus 
tidal perturbations. This latter approach has proved to be very powerful and we have succeeded in obtaining the 
same equation of motion as obtained in the former approach. 

The correction term of 0{Gm) is found to be entirely given by the part of the metric perturbation which is due 
to the tail term of the Green function. Defining the regularized metric as the background plus this tail part of the 
perturbation, we have found that the equation of motion is the geodesic equation on this regularized perturbed metric. 

It is important to note that the particle does not have to be a black hole but the resultant equation of motion can 
be equally applicable to any compact bodies such as neutron stars. The essential assumption here is that the only 
scale associated with the particle is Gm and the structure is basically spherically symmetric. In this sense, we have 
shown the strong equivalence principle to the first order in Gm. 

We also emphasize that our result gives a justification of the so-called black hole perturbation approach for the first 
time. In the black hole perturbation approach, one calculates the gravitational radiation from a particle orbiting a 
black hole with the assumption that the particle is a point-like object with the energy momentum tensor described 
by the delta function. Although this approach has been fruitful, there has been always skepticism about the validity 
of the delta functional source. What we have shown in this paper justifies the use of the delta function in the source 
energy momentum tensor. 

In this paper, we has considered a particle which is essentially structureless. A next step will be to include the 
intrinsic angular momentum of the particle, i.e, to consider a spinning particle. As we have found the method of the 
matched asymptotic expansion is very powerful, there is a big possibility that we can extend the present analysis to 
the case of a spinning particle. Extension in this direction is under study. 



We thank T. Futamase, T. Nakamura, M. Shibata, H. Tagoshi, K.S. Thorne and A.G. Wiseman for fruitful con- 
versations. Y. M. thanks Prof. H. Sato and Prof. S. Ikeuchi for their continuous encouragements. This work was 
supported in part by Monbusho Grant-in- Aid for Scientific Research No. 5427. 



Here in this appendix, we summarize our notation and basic formulas. As for the formulas already derived in DB, 
we just write down the results. 



VI. CONCLUSION 
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APPENDIX A: NOTATION AND BASIC FORMULAS 



1. Basic Notation 



(i) The Riemann tensor and the Ricci tensor are defined by 




+ r^^.^p'^ 




(Al) 
(A2) 
(A3) 
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(ii) Symmetrization and anti-symmetrization of the tensor indices are described by t^^") = {tf" + r^)/2 and t^f"^ = 
^-l-f,^ _ t''f^)/2, respectively. 

(iii) For an arbitrary bi-tensor Q, 

Q;a{x,z{Ta;)) := [Q;^ (x, , 

Q;ti {x,z{Ta;)) [Q ■f.ix , z)]^^^(r^) , (A4) 

while 

[Q{x, z(Ta;))];^ Q-.fj^ix, z{tx)) + Q-a{x, z{t^))z°'{t^)tx.^. (A5) 

(iv) The basic equations satisfied by a half the squared geodetic interval a{x, z) (Eqs. (1.11) and (1.12) of DB): 

(j{x,z) = ^g^"" {x)a.^{x, z)(T.„{x, z) = ^g°''^{z)a.a{x, z)(T.j3{x, z), 

lim a{x, z) = 0. (A6) 

X — 

(v) The defining equations of the parallel displacement bi-vector (Eqs. (1.31) and (1-32) of DB): 

gtia;i^{x, z)g'"^{x)a.cr{x, z) = 0, gi_ta:f3{x, z)g^^ {z)a.^{x, z) = 0, 

lim5^"(x,z) = V- (A7) 

(vi) The definition of A(a;,z) (Eqs. (1.50), (1.51), (1.60) and (1.61) of DB): 

A(x,z) ■.^\r^{z,x)a,,,0{x,z)\. (A8) 



2. Basic formulas 



(i) (Eq. (1.28) of DB) 



(ii) 



a-afiix, z) = gafiiz) + '^RoP p^[z)<7.^{x, z)<7.5{x, z) + O(e^). 



r'^A^^ ^) = \r''{x, z)R",0s{zy''{x, z) + 0{e^), 
r'^A^, z) = \r\x, z)g,^x, z)R^p^s{z)a-\x, z) + 0{e^), 



which can be obtained from Eqs. (1.40) and (1.41) of DB. 
(iii) (Eqs. (1.51), (1.64) and (1.73) of DB) 



(A9) 



(AlO) 



(J-,j.p{x, z) = -gi_"{x, z) ( gafiiz) - -Ra-,i38{z)(j'''^ [x, z)a- (x, z) ) + 0(e ). 



(All) 



(iv) Here we give the formulas which we n eed t o expand Eq. ( ^.20 ) to obtain Eq. (2.27). Since we want to express 
Eq. (2.20) in terms of defined by E q. (^.23 ) instead of Tjfet/Advix), we expand each factor of each term which 
consists of Eq. ( 2.20 ) with Sj^et/Advix), ( ^.26| ). We first consider ['^{x, z{T))]^^^^^^^^^^f^^s^, which is expanded as 



[cr{x,z{T))]^__ 



■T Het/Ad 



„(x) = ^i^' ^i'^x)) + (^{x, z{T^))5Ret/Adv{x) 



1 .... 



+ - a {x,z{T^))Sji^t/A^^{x) + — "& (a;,z(T^))(5fl^t/^rf„(x) + 0(e ) 



(A12) 



Each term is computed as follows; 
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= a-afsix, z{Tx))z°'{Ta,)z'^{Tx) + a-aix , z{Tx)) z" (t^) 

gap{z(Tx)) + ^Ra'' p\z{T^))a.^{x,z{T^))a.s{x,z{T^))^ z" (t^) z'^ (t^) 

+a,Jx, z{t^))z°'{t^) + 0{e^), 

a (x,z{t^)) = a.a{x,z{T^)) "z °'{t^) + O(e^), 
a {x,z{t^)) = -5Q/3(^(Tx))^"(T:r)i;^(ra;) + 0(e). 



(A13) 



In the above computation, we used Eqs. (A9), ( 2.23 ), and the normalization condition (dz/dr)^ = —1 
which was proved to be consistent in sections 3 and 4. Then we obtain 



1 



ct(x,z(t)) 



± 



T = Tjiet/Atii,(2;) 



^^^1^ (^1 T le{x) z "iT,)a,^{x, z{t.,)) - le\xyP{T.,) + 0{^) 



Noting the exphcit form of u^^°'^{x,z) in Eq. ( 2.1C ), it is necessary to compute [A^/^(a;, z(t))]t- 
[gMa(^'^W)]r=r„et/Ad„(^) and [i^"(T)]r=r„,t/^^,(x)- In the same way, 



(A14) 
(A15) 
(A16) 

0(Gm/i), 



(A17) 



= Tj?.et//l£ii,(3:)> 



^^/\x,z{t)) 

gf,a{x,z{T)) 



= l + 0(e'), 



T=THct/Adi,(a:) 



■r=Tiict/Adv(x) 

(Jf^aix, z{t^)) ± ^g^'^(a;, z{Ta:))Ra0^s{z{Tx))a''^ {x, z{Ta;))z\Tx)e{x) + O(e^), 



i"(r) 



= i"(T,) T e(a;)K-i(x)z"(T,) + -e\x) z "(r,) + Oie^), 



Putting them into the first term in the parentheses of Eq. (2.2C), and using Eq. (2.11) and the fact 

A'^ix) = ro.{x, z) (^^"(z) - a,p{x, z)A^-^\z) + ©(e^)^ , 



(A18) 

(A19) 
(A20) 

(A21) 



in computing the second term, we obtain Eq. (2.27) 



APPENDIX B: 2ND ORDER VARIATION OF THE EINSTEIN TENSOR 



We derive Eq. (3^) by taking the variation of the Einstein- Hilbert action. We first compute the Einstein- Hilbert 
action of the metric, g^j^y — g^^ -V h^^. For this purpose we define a differential operator Sg: 



SgQ{g) = lini 



Qig + eh)-Qig) 



e^O e 

.2n( ^ r Q{g + ^h) + Q{g - eh) - 2Q{g) 
5gQ{g) = lim . 

We first note the 1st variation of the Einstein tensor taken from a standard text book. 



(Bl) 
(B2) 



(B3) 
(B4) 



where -tpp^, = hp^^, - {l/2)gpyh. 

Since we are interested only in the 2nd order variation, we compute the terms proportional to in the action: 
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where we have used the assumption that the background spacetime is vacuum, i?^^ — 0. Now 



Inserting these into Eq. (B5) and using 



1 



we obtain the 2nd order variation of the Einstein-Hilbert action, 

1 1 



(B5) 



(B6) 
(B7) 



(B8) 



-SliV^R) 
IdttG V — (? 



1 



647rG 



where ^p = ip^a = — ^ and we have discarded unimportant total divergence terms. 
We note 



SS SS , 

J^fe] = T I 



(B9) 



(BIO) 



Thus the 2nd order variation of the Einstein tensor can be obtained by taking the variation of the action with respect 
to g^u- 



Hence, 



ogp,^ 



(Bll) 



(B12) 



Carrying out the variation of L^^^ one finds 

'1 



1 



(B13) 



Note that if h is a hnear perturbation which satisfies G^^'>'"'[h] = 0, the tensor T^" -(l/87rG)G(2''''' describes the 
conserved energy momentum tensor of the perturbed gravitational field. 
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APPENDIX C: TENSOR HARMONICS EXPANSION 



Here we briefly review the construction of tlie scalar and the vector harmonics in terms of the symmetric trace-free 
(STF) tensor We introduce the notation 



A 



<iii2---ip_> 1 



(CI) 



to represent the totally symmetric and trace-free part of Aiji^...^^. More explicitly in the cases of ^ = 2, 3, 
The spherical harmonics expansion of a scalar function A on the unit-sphere can be written as 

oo 

A = A^i-^i 



(C2) 



(C3) 



£=0 



where — X'^/\X\. In this case, the order £, which is associated with the angular dependence, is equivalent to the 
total angular momentum, J. Thus the J mode of the (TT)-component of the metric perturbation is totally determined 
by its angular dependence. Namely, the terms in the (rr)-component of the metric perturbation which contain 



1, n\ n^'n^^, 



(C4) 



correspond to the J = 0, 1, 2 modes, respectively. 
Next we consider the expansion of a vector field Ai , 



(C5) 



In this case the term of the £-th order in the angular dependence is decomposed into J = £ + 1, I and This is 

done by using the Clebsch-Gordan reduction formula |l4| , 



= R 



(+) 



R 



(0) 



2^-1 



Si^i.R. 



(-) 



(C6) 



where Ti^i^...i^ is a STF tensor of order £ and 



(+) 



— C^<'if+i?iii2 ■■■»{>) 



(C7) 



We perform the decomposition explicitly for £ < 2 here. For £ = 0, there exists no J = mode and it trivially 
corresponds to the J = 1 mode. For £ = 1, the decomposition is performed as 



Ai 



ls,,A^ 



ij^kk 



A[ij] + -SijAkk 



(C8) 



and the first, second and third terms in the square brackets correspond to the J = 2 ,1 and modes, respectively. 
For £ — 2, we obtain the decomposition formula as 



A,<jk>n<'n''> = 



A 



(2) 



(C9) 



where 



(2) 



(CIO) 
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and the first, second and third terms correspond to the J = 3, 2 and 1 modes, respectively. 
As an example, let us consider a vector, 

dz°' 



(Cll) 



which appears in the (q) and (^) matchings in section 4. First we decompose it in terms of its angular dependence as 

(C12) 



Using the relation ( 4.16| ) and the fact that the Ricci tensor vanishes, the second term in the parentheses is rewritten 

1 dz°' „ „^ dzf^ dz^ , 



as 



3 dT 



■R^0,sf\^^\X\ 



(C13) 



and is found to be zero due to the symmetry of the Riemann tensor. Thus we have onlv to c onside r th e first term in 
the parenthesis of Eq. (C12), which is decomposed further with the aid of the formulas ( |C10 ) and (C£) as 



dz" 

It 



(C14) 



where 



(2)7/35 _ 



r, {,f^mf <nf i>^jmn ~t~ f^mf <nf 



(1)7,35 



(C15) 



It is easy to see that the first and third terms vanish due to the symmetry of the Riemann tensor and the Ricci 
flatness. Thus only the J = 2 mode remains. 

As for a tensor field, it is not necessary for us to give general discussion here. The only term which requires our 
consideration is 



D ^ct fB /•7 fS \rK \^ra 

n-ap-iSj ij kj J J ^ , 



(C16) 



which appears in the (zj)-component of the (q) matching. In this case, it is better to use the symmetry of the Riemann 
tensor from the beginning. First, we define the spatial triad components of the Riemann tensor by 



Rijkm ■ — Rap-ySf if^jf^kf rn ■ 

Introducing a symmetric tensor defined by 



we can express R'^^^'^^ in terms of TZij as 



~~ ^ ^ikm^jns Rkn 



p nij skm^ 



(C17) 
(C18) 

(C19) 



Then the symmetric tensor TZij is decomposed into STF tensors as 



TZij — 7?.<y> + —SijTZkk- 



(C20) 

Counting the number of indices, we find that the first and second terms in Eq. (C2C) correspond to J = 2 and 
modes, respectively. However, again owing to the symmetry of the Riemann tensor and the Ricci flatness, the J — 
mode vanishes and only the J = 2 mode remains. 
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